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Abstract 

We consider a Hofava theory that has a consistent structure of constraints and 
propagates two physical degrees of freedom. The Lagrangian includes the terms 
of Bias, Pujolas ans Sibiryakov. The theory can be obtained from the general 
Horava's formulation by setting A = 1/3. This value of A is protected in the 
quantum formulation of the theory by the presence of a constraint. The theory has 
two second-class constraints that are absent for other values of A. They remove the 
extra scalar mode. There is no strong-coupling problem in this theory since there 
is no extra mode. We perform explicit computations on a model that put together 
a z = 1 term and the IR effective action. We also show that the lowest-order 
perturbative version of the IR effective theory has a dynamics identical to the one 
of linearized general relativity. Therefore, this theory is smoothly recovered at the 
deepest IR without discontinuities in the physical degrees of freedom. 



1 j orgebellorinOusb . ve 
2 arestu@usb . ve 
3 asotomayor@uantof . cl 



1 Introduction 



The consistency of Hofava theory [T] is a subject that has been under intense study 
motivated by the search of a perturbatively renormalizable theory of quantum gravity. 
Numerous models that follow the original proposal of Horava of adopting the foliation- 
preserving diffeomorphisms as gauge symmetry have been analyzed. An unavoidable 
question for all these models is whether general relativity (GR) can be consistently re- 
covered at large distances, such that the Horava theory can be regarded as the UV com- 
pletion of GR. The originally proposed scheme for emergent GR inside Horava theory [T] 
consists of obtaining the full action of GR (in the Arnowitt-Deser-Misner formulation) 
as the lowest-order effective IR action. This scheme requires the enhancing of the gauge 
symmetry group to the general space-time diffeomorphisms; that is, the original proposal 
is that not only the dynamics but also the symmetry of GR should be restored from 
Horava theory in an approximate, effective way. To achieve this, the IR effective action 
should satisfy the condition that the coupling constant A multiplying the trace-kinetic 
term, 



must approach the value A = 1 in order to this combination matches its fully covariant 
version. 

An alternative to this scheme was noticed in Ref. [2], where it was shown, by means of 
a Hamiltonian analysis, that the lowest-order truncation of the original, nonprojectable, 
Horava theory, which is given by the Lagrangian 



is physically equivalent to a gauge-fixed version of general relativity (the gauge in which 
K = 0, the so-called maximal slicing gauge). This means that both theories are dynam- 
ically identical, although their gauge symmetry groups are different. The crucial point 
for the result of Ref. [2] is that the condition K = emerges as one of the constraints of 
the theory, hence A becomes meaningless for the IR effective action. In particular, this 
result shows how it is possible to get GR without the requisite A — > 1. 

However, there is a central question about the discontinuities that might arise in 
recovering GR. Since Hofava theory has a reduced gauge symmetry group, generically it 
has an extra degree of freedom with respect to GR (this was already studied in Ref. [1]). 
There is an abundant quantity of works devoted to the physics of this extra mode. We 
may mention that in Ref. [3] it was signaled the problem of its strong coupling in the 
original theory; implying, instead of a discontinuity, the breaking down of the whole 
perturbative analysis. We want to stress that this result is based upon the assumption 
that A — > 1 at the IR as a generic rule to get GR. In Ref. [3], using a curvature-square 
model, it was shown that the extra mode is of odd nature, that is, propagates itself with a 
first-order time derivative. This was confirmed in Ref. [6], using again a curvature-square 
model, but showing also that the algebra of constraints closes and that the extra mode 
decouples at large distances smoothly in perturbative analysis^] 0. 

4 The results of Refs. j2J[5] were corroborated in Ref. [TJ. 

5 A perturbative analysis in a projectable model of the theory can be found in Ref. [5]. 
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With the aim of curing the oddness of the extra mode, it was noticed in Ref. [9] that 
the nonprojectable Hofava action admits a large class of new terms (once the principle 
of detailed balance is discarded). Since the vector a, = dihaN is covariant under the 
foliation-preserving diffeomorphisms, scalar combinations of it and the spatial metric 
are admissible into the Lagrangian. Adopting the logic of renormalizable gauge field 
theories, all the terms that are compatible with the gauge symmetry must be included in 
the Lagrangian. This leads us to consider the complete, nonprojectable, Hofava theory as 
the one containing the terms of Bias, Pujolas and Sibiryakov [9|. Those authors found in 
[9] that the extra mode becomes even (propagates with a second-order time derivative) in 
the complete theory. However, the authors of Ref. [10] reported that the strong coupling 
problem persists in the complete theory, assuming again the condition A — > 1 at the IR 
(see also [H]). It has been argued p2] that this problem can be avoided by requiring 
that the scale at which high-order operators become relevant is low enough. 

On the side of the Hamiltonian analyses, in Refs. [131 HH EH] it was studied the closure 
of the algebra of constraints of the complete theory. In Refs. [131 [Tj)] explicit computations 
were performed on the lowest-order IR effective action (second-order in derivatives). It 
was shown that the algebra of constraints closes: in particular well-behaved elliptic partial 
differential equations arise (see [IS]), and these analyses confirm that the theory has an 
extra even scalar mode, which is persistent at the level of the IR effective action. 

It is important to realize that there is an underlying hypothesis behind these Hamil- 
tonian analyses: it is assumed that the time derivatives g^ can be completely solved in 
terms of the momenta 7r u [15]. This is effectively the case when A satisfies A 7^ 1/3. 
However, when A = 1/3 it is not possible to solve completely in terms of 7r* J ; instead, 
the primary constraint g^n^ = emerges. Since this constraint is not present in the 
case A 7^ 1/3, it is to be expected that the number of physical degrees of freedom is 
reduced when A = 1/3. Therefore, it becomes of great interest to study the Hamiltonian 
formulation for the complete Hofava theory under the special value A = 1/3, since this 
case might be an exception to the generic presence of the extra mo Moreover, one 
may ask whether the dynamics of theory at A = 1/3 is closer to the one of GR than 
for other values of A. We emphasize that the gijir^ = constraint protects the value 
A = 1/3, since any other value for A would imply the violation of this constraint and the 
quantum formulation must be done on the constrained submanifold. 

With these goals in mind, in this paper we perform the Hamiltonian analysis to 
the complete, nonprojectable Horava theory fixing the special value A = 1/3. Explicit 
computations on the complete Hofava theory are very difficult since the Lagrangian has a 
big number of higher-order terms. To overcome this difficulty, our first strategy consists of 
dealing with a general potential and making computations in an implicit form. This will 
allows us to arrive at conclusive results on the closure of the algebra of constraints and 
the number of physical degrees of freedom. Then we move to a specific model which has a 
z = 3 term and the most general z = 1 terms. The z = 3 term is the square-Cotton tensor 
term and can be obtained by the detailed balance principle [U [IB]- The z = 1 terms give 
the relevant action for the large-distance dynamics, the IR effective action. Of course, this 
is justified by the assumption that all higher-order terms are suppressed at low energies. 



6 The Hamiltonian formulation of the lowest-order truncation of the original, nonprojectable, Hofava 
theory at the value A = 1/3 was carried out in Refs. p~7| 12]. 
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This model, which will allow us to perform explicit computations, can be regarded as 
a theory with soft breaking of conformal symmetry, because the square-Cotton term is 
conformally invariant whereas the z = 1 terms break the conformal symmetry. 

We give in advance our three main results: The Hofava theory we consider has a closed 
structure of constraints, propagates two physical degrees of freedom and the linearized 
version of its IR effective action coincides with linearized GR. Notice that these results 
imply that linearized GR is recovered at the lowest energies without discontinuities in 
the physical degrees of freedom. For the concrete model we analyze we also found that 
in a sector of the space of parameters the energy of the model is nonnegativeQ 



2 Hamiltonian analysis of the full theory 

We now start with the computations. Most of the steps are parallel to the ones of the 
case A 7^ 1/3. In order to have a self-contained study, we shall perform the whole analysis 
from the very beginning, making special emphasis on the results that depart from the 
A 7^ 1/3 case. 

The action of the complete, nonprojectable Hofava theory is written in terms of the 
Arnowitt-Deser-Misner (ADM) variables gij, N and N{ as 

S = J dtd 3 x^N {G'^K.jK,, - V) , (2.1) 

where 

K* = ^(^-2V ( ^)), (2.2) 

G ijkl = \(g ik g jl + g a g jk )-Xg ij g kl (2.3) 

and the potential V = V(gij, o«, • • •) is the most general combination of the spatial met- 
ric, its curvature tensor, the vector Oj and covariant spatial derivatives of these objects 
that transforms as a scalar under spatial diffeomorphisms. To ensure power-counting 
renormalizability, the potential must include at least terms of order z — 3, which means 
that they are of sixth order in spatial derivatives. The potential can also include a 
cosmological-constant term; we put it equal to zero, hence our simplest vacuum is the 
Minkowski space-time. The lowest-order terms, which yield the effective action for large 
distances, are 

V (2) = -R - aaia* . (2.4) 

A and a are coupling constants of the theory. 

By expanding the kinetic term we get G i3 K^K^i = K^K 1 ^ — XK 2 , where K = g^R^. 
If we were dealing with GR we were forced to put A = 1, since only the combination 
KijK l i — K 2 is covariant under transformations mixing time with space. However, under 
foliation-preserving diffeomorphisms, both K^K 1 ^ and K are separately covariant, hence 

7 Positiveness theorems for the energy of nonprojectable Hofava theory in the case A ^ 1/3 have been 
formulated in Refs. p~6 | 1X9 ] . 
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A is left undetermined by the gauge symmetry of the Hofava theory. From the general 
relation 

G^ gk i = (l-3\)g^, (2.5) 

we may see that the value A = 1/3 is special since the metric becomes a null eigenvector 
of the four-index metric G u7fc , 

G^ lk g kl = 0. (2.6) 

This implies that the metric G l ^ kl is not invertible for A = 1/3. 

Let us perform the Legendre transformation to the action (12.1 1) in the case A = 1/3. 
Since general spatial diffeomorphisms are part of the gauge symmetries of the theory, we 
know that the shift functions iVj can be regarded as the Lagrange multipliers associated 
to the first-class constraint generating this kind of transformations. Hence the phase 
space is spanned by the conjugated pairs (g^, 7r 4J ) and (N,<p). The action (12 .ip does not 
depend explicitly on the time derivative of N, hence we get the primary constraint 

= 0. (2.7) 

The momentum conjugated to the spatial metric is given by 

^l = G^ kl K kl . (2.8) 

By using ( 12. 6ft and denoting tt = gijit l \ we get the primary constraint 

tt = 0. (2.9) 

This constraint is absent in the complete, nonprojectable Horava theory with A ^ 1/3 
[131 [TU [15]. 7r is the generator of conformal transformations on g^ and ir^ . However, 
these are not part of the gauge symmetries of the theory, hence we may anticipate that 
7r = is a second-class constrainto A similar consideration applies for = 0. 
From ( 12. 8ft it is straightforward to obtain the relations 

G ljkl K tJ K kl = -tt^, (2.10) 

= ^~ U ~<j ■ - , ~' / V,.V / . (2.11) 

and using them we may build the Hamiltonian. We get, after an integration by parts 
that yields no boundary contributions, 

H = J d 3 x (^K ij TTij + v^iVV + NiW + a<P + fii^J + E ADM , (2.12) 

where "H* = is a primary constraint defined by 

T-C = -2V j vr ij + 0N . (2.13) 



8 For the special case of a potential containing only the square of the Cotton tensor the action acquires 
a conformal symmetry at A = 1/3 [1] . 
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This is the so-called momentum constraint, which is the generator of spatial diffeomor- 
phisms, hence it is a first-class constraint. The term proportional to (j) in T-L l ensures 
to have the complete generator, since in the complete Horava theory N is part of the 
canonical variables. We may add this term since <fi = is a constraint of the theory. We 
have also added the rest of primary constraints to the Hamiltonian (12.121) . such that iVj, 
a and /i enter as Lagrange multipliers. Obviously, fi does not arise in the case A ^ 1/3. 
The ADM energy, 

-Eadm = j> d£i{djgij - digjj) , (2.14) 

is included a la Regge and Teitelboim [20] in order to obtain the equations of motion from 
the most general variations of the Hamiltonian that are compatible with the boundary 
conditions. 

Now we study the preservation in time of the primary constraints. Since W = 
is a first-class constraint, we concentrate ourselves in the preservation of <fi = and 
7i = 0. For the time evolution of <ft we need to compute its Poisson bracket with the 
Hamiltonian (12.121) . In computing the Poisson brackets we can omit the boundary terms 
of the Hamiltonian, since they do not give local contributions to these brackets. We just 
need to remember that we may discard any nonzero boundary contribution arising in the 
derivative of the bulk terms, since the boundary terms of the Hamiltonian account for 
them. Notice that arbitrary variations of the potential with respect to N can be written 
in a closed form, 

S N J d 3 x^NV = j d 3 x^/g (vSN + N^S ai + N d( ^. a ^ J a i + ' ' " J ( 2 - 15 ) 

= [d 3 x^l V5iV + iV^— s V ir ... i2 8a h ) , (2-16) 

J y r=1 0{V ir ...i 2 a il ) j 

where we are using the shorthand V^ -. = VjVjVfe • ■ •. After using 5a« = di(SN/N) and 
integrating by parts (no term gives boundary contributions), we get 



j d 3 x^NV= I d 3 x^5N y+-L]T(-irV n ... v ( 

J J r=l V 



N- dV 



(2.17) 



d{V ir ... i2 a h 

Therefore, the condition {<p, H} = leads to the secondary constraint Ti — 0, where 

K^-Ltt^. + ^V (2.18) 
and we have introduced the modified potential 

* = V + 7j E(- 1 ) rv - - ( N f>m dV s ) ■ ( 2 - 19 ) 
Following the standard nomenclature of GR, we call % = the Hamiltonian constraint. 
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The Hamiltonian given in f 1 2 . 1 2 j) can be written in terms of the other constraints 
and boundary terms. For, we notice that an integration by parts and the behavior of the 
fields at infinity allow us to verify the identity 

J d 3 xNU = j d 3 x (^Tr^TTij + y/gNV^j + 2a$ N , (2.20) 

where <3>tv is the flux of N at spatial infinity, 

$tv = {CLEAN. (2.21) 



To arrive at the integral (I2.20p we have used the fact that all the derivatives of the 
potential enter in N% inside total divergences; the only one that does not vanish when 
integrated at spatial infinity is the derivative of the quadratic term in Oj, which yields 
2a^fgNd l = 0(r~ 2 ). Its surface integral at infinity can be further simplified such that 
it gives the flux of N. On the basis of (I2.20p we may write the Hamiltonian ( I2.12p as a 
sum of constraints plus boundary terms, 



H = 



J d 3 x (NU + NiW + (70 + fm) + E ADU - 2a$ N . (2.22) 

From this expression it is clear that the energy of the theory is given by 

E = E ADM - 2a$ N . (2.23) 

From a mathematical point of view, the role of the boundary term proportional to $jv 
in the Hamiltonian f 1 2 . 2 2 j) is analogous to the one of the ADM energy (20] : it ensures the 
differentiability of the Hamiltonian under variations of N that behave as SN = 0(r~ x ) at 
infinity. In particular, we shall see below that there is a term proportional to V,a* in the 
modified potential V. Its variation with respect to N gives rise to a nonzero boundary 
term that cancels out with the variation of — 2a$7v- 

We now turn our attention to the preservation of n = 0. We find that the condition 
{%, H} = leads to the constraint C — 0, where 

3iV 

C=— tt^-V^V', (2.24) 

and we have introduced the objects 

Hj = ^J d 3 y^NV , V ' = 9ij V Hj . (2.25) 

Again, this constraint is absent in the complete theory with A ^ 1/3. Notice that there 
are some similarities in the structures of H and C. Two immediate consequences we 
obtain from the vanishing of them are 

V<0, V' = --NV. (2.26) 
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The next step is to impose the preservation in time of the secondary constraints 
% = and C = 0. By computing Poisson brackets with the Hamiltonian we obtain that 
the condition {H, H} = leads to the equation 

/ *y *± J d 3 z^V5 + J d *y, 9lj A / d 3 z^V5 - ^ ^ 

+2 f d 3 y^^- f d 3 z^gVS-2^V' ij = 0. 

The symbol 5 is the Dirac delta centered at the point x at which this equation is evaluated, 
5 = 5 3 (z — x) 0. Similarly, the condition {C, if} = leads us to the equation 

J y/g Sgtj J 

(2.28) 

Eqs. (I2.27ti2.28j) form a system of equations for the Lagrange multipliers a and fi. 
Indeed, since the potential depends on derivatives of iV and <7y, the first two terms of 
both (I2.27P and (I2.28|) lead to differential operators on these multipliers. Thus, Eqs. 
(I2.27H2.28"]) are a coupled system of partial differential equations (PDEs) for a and /i. 
Whenever this system can be solved for a and /i, Dirac's algorithm for the preservation 
of constraints ends consistently with these equations (in the complete theory with A 7^ 1/3 
Dirac's procedure ends with a PDE for a [HI IT5]). 

Moreover, for the consistency of the whole Hamiltonian formulation it is of central 
importance that the Eqs. f)2.27H2.28l) can be solved for the multipliers a and /1, without 
any further restriction on the canonical variables. Since the most relevant terms to 
determine the existence of solutions are the highest-derivative terms, it is illustrative to 
study the structure of the highest-order terms that arise in Eqs. f)2.27H2.28l) when the full 
(up to z — 3) potential of the theory is considered. There is a big number of inequivalent 
z = 3 operators that can be constructed with the curvature tensor and af, thus, at first 
sight it seems a very difficult task to elucidate the structure of the highest-order terms of 
Eqs. fl2.27H2.2Hl) . However, a direct analysis on these equations may convince ourselves 
that some of these operators lead to the cube of the Laplacian, V 6 , acting on a or /i and 
the other ones lead to lower-order terms. For, the main point is that in Eqs. f)2.27H2.28l) 
the four terms 



J d3ya lN /^^' / ' d3 y^9ij^~ J d 3 Zy ^V5, 
J d 3 y^g tJ ^- J d 3 z^V'5, J d 3 ya-^ J d 3 z^V'5 



(2.29) 



Let us further explain the notation. For example, the first term of (|2 . 2T[) should be read as 



5N(y) 

and similarly the other terms of (|2.27|) and (|2.28|) that involve integrals and functional derivatives. 
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contain derivatives of the potential V of at least second order and there are z = 3 
operators that combine all their derivatives into a sixth-order derivative when they are 
functionally derivated twice (or more). Let us consider, for example, the two operators 
£>i = {ViR jk ) 2 and 2 = V i a i V 2 V j a j . When O x is inserted in the third term of $£M), 
the highest derivative that it gives rise on the Dirac delta is V 6 . After integration by 
parts, we get that 0\ yields a term proportional to V 6 /i in Eq. (I2.28p . whereas all its 
other contributions are of lower order. Similarly, O2 yields V 6 acting on the Dirac delta 
when inserted in the first term of (12.291) . hence it leads to a term proportional to V 6 <7 in 
Eq. ( 12.271) . The operators 0\ and O2 must be included in the potential in their direct 
forms or in terms of their equivalent operators; that is, other operators that are obtained 
from them by integration by parts. In any case we obtain V 6 as the highest-order operator 
acting on the Dirac delta once the functional derivatives in (12.291) are performed. Other 
inequivalent operators yield the same result when inserted in some of the terms in (I2.29P . 
We could anticipate this result by noting that there is no covariant differential operator of 
sixth order acting on the Dirac delta other than the cube of the Laplacian. Note, however, 
that there are z = 3 terms that do not yield sixth-order derivatives acting on the delta 
when inserted in (12.291) . but lower-order derivatives. An example is (a^a*) 3 . Another 
interesting example is the square Cotton, C^Cij, which we shall consider explicitly in the 
next section. 

In despite of the fact that some operators yield lower-order derivatives on a and fi in 
Eqs. (I2.27H2.2H1) . the ones that yield V 6 must be included in the potential and these are 
the dominant ones in Eqs. (I2.27tl2~28|) . After all the terms proportional to V 6 a and V 6 // 
are collected, the last requisite is to impose that the matrix of their coefficients is positive 
definite, which is a condition in the space of coupling constant^. Thus, we have that, 
when the full potential of the theory is considered, Eqs. (l2.27tE2T2"%j) constitute a system 
of sixth-order, linear, elliptic PDEs for a and /x characterized by the highest-order terms 
VV and VV 

Having seen that the theory has a closed structure of constraints, we now evaluate the 
number of degrees of freedom. The theory has the first-class constraint W = and the 
second-class ones <p = 7i = 'H = C = 0. They leave four independent degrees of freedom 
in the phase space, which correspond to the propagation of two even physical degrees 
of freedom. This is the same number of GR; there are not extra degrees of freedom in 
the complete, nonprojectable Hofava theory at the value A = 1/3. We may regard the 
Hamiltonian constraint "H = and the constraint C = 0, where % is given in ( 12 . 1 8[) and 
C in (I2.24p . as a coupled system of PDEs for the lapse function N and one mode coming 
from the spatial metric g,ij. 

10 Notice that there is not place for nonconstant coefficients in these terms. The potential only depends 
explicitly on gij and at and their derivatives; any coefficient that depends on them necessarily would 
increase the order, and we are considering the highest-order terms. 
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3 Soft breaking of the conformal symmetry 
3.1 The full model 

In order to develop explicit computations, in this section we consider a model with z = 3 
and z = 1 terms. This model will help us in clarifying the structure of the above equations. 
The quadratic (z = 1) terms we consider are the most general ones; they are grouped in 
V (2) , which is given in ( 12. 4p . The importance of this quadratic potential lies on the fact 
that it is the leading one at the lowest energies, such that we build the lowest-order IR 
effective action with it and the kinetic term. This is the appropriated scenario to test 
whether GR is recovered at the IR. For the z = 3 term we consider the square-Cotton 
term, which was the original z = 3 term proposed by Hofava [T] and can be elegantly 
justified in the 3+1 action by the detailed balance principle. 
The potential is given by 

V = -R- aa ia l + «•<",/ "' , (3.1) 
where w is a coupling constant and C iJ is the Cotton tensor, 

& = — e kl ^V k Ri j) . (3.2) 

Computations with the C 2 term are facilitated by the fact that JgC^Cij transforms 
homogeneously (with weight —3/2) under conformal transformations of the metric, which 
are generated by it. 

The modified potential (I2.19P and its derivative V ( I2.25P take the form (before im- 
posing constraints) 

V = -R + a (2V i a i + a^) + wCijC ij , (3.3) 

V' = --NR + 2NV i a l + (2-a/2)Na i a i -^NC ij C ij , (3.4) 

such that the Hamiltonian constraint (12.181) and the C constraint (12.241) become 
U = -^~' J ~,j - y/gR + a^/g (2V l a i + a i( j) + w^gC % f ij = , (3.5) 

C = ^^7r t , + ^NR-^N(2V l a t + (2-a/2)a l a l ) + ^^NC lJ C i ^ = 0. 

(3.6) 

Finally, we evaluate the equations (12. 2 7112.28]) for the model ( 13. II) . As we anticipate in 
the previous section, these equations lead to a coupled system of PDEs for a and [i, 

f3 (2VV + Ntfdiii) - 2g- 1 ii ij ix ij a + (f3V 2 N + 3// ; .V;t' ; ;,•,,) /i + (3^A^ - 2a)wC ij C ij = 

AN 4/3 
irVfNRa - ViVjN + aNaidi) + —dANd^N^) 



wN „ s 8wN 

ij) i 



djO^iN^) + -—<K kl &> kl {NC K 



V9 J V9 
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(3.7) 



V 2 yU - ^a l di<7 - - (R + aaia 1 + (3/j)g l n lj 'K ij ) fi + ^a { a l a - ^-CijC lj \x 
2a ,, „ „ 2w „ v , 2w _.. w 



-tt ij (iVi? - ViVjiV + aiVa^) C O^A^) + -— 7r fc/ CF fc <(^) , 

(3i 

where O lJ ' fcJ and 0^ fe are the differential operators 
m _ 1 



_ (^ mfc V m V V + g jl e imn V m V k V n - g^e imk V m V n V n - g kl e imn V m V j V n ) 

+-^= (g jk e iln R n m V m - g^ k e imn R n l V m - g^e lkn R n l V m ) , 

W kl = -L {-e imk V j VV m - g jl e imn V n V k V m + g jl e imk V n V n V m + g kl e irnn V n V j V m ) 
2^g 

+^V m [(-g jk e Un R n m + gi k e %mn R n l + g^ m e ikn R n l ) ■] , 

2\J g 



and 



(3.9) 

^(r^)' (310) 

We have made some massaging on Eqs. 03.7H3.8p to brought them to their final form, in 
particular by using the constraints (I3.5H3.6I) . Eqs. (l3.7H3.8p constitute a system of linear, 
elliptic, PDEs for the Lagrange multipliers a and fi. Notice that these equations are 
of second order although the model is z = 3. This is because the JgNCijC^ term is 
covariant under conformal transformations generated by N<fr and n. Thus, the Poisson 
brackets between yfgNCijC 1 ^ and acj), fiir generate terms proportional to the C 2 term, a 
and \x playing the role of conformal factors. Hence no derivatives of a or /i are generated 
by the C 2 term. Eqs. ( l3.7H3.8p can be solved for a and \x on general grounds, hence we 
conclude that Dirac's algorithm for the preservation of constraints ends consistently with 
these equations. 

As we pointed out in the previous section, constraints (l3.5H3.6p can be regarded as a 
system of PDEs for A and one mode coming from gij. We would like to emphasize this 
point for N. Constraints (l3.5H3.6p can be combined to obtain the equation 

(3S/ 2 N = A Q^TTy + wC ij C ij J (3.11) 

The structure of this equation falls on the same class of equations we studied in Ref . [16] . 
There we showed that, if the term without derivatives of A satisfies a positiveness condi- 
tion, then the solution for N, in the sense of distributions, exists and is unique. For the 
Eq. (13. lip , the nonderivative term is definite positive by requiring only w > 0. Assuming 
also a < 2 (j3 > 0), then the results of Ref. [T6] imply the existence and uniqueness of A 
and moreover it is guarantied that A > over all the spatial submanifold. 
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By evaluating conditions (I2.26P (or directly from H3.5ti3.6j) ) . we get, imposing w > 0, 



R- (l + 3a/2)vV - (l + a/2)a i a i = 0, (3.12) 
R - 2aV i a i - aa^ > 0. (3.13) 

It is straightforward to deduce from these relations the following inequalities 

P V 2 iV > , 7 (R + acne?) > . (3. 14) 

On the basis of the inequalities (I3.14p . we may give a result about the positiveness of 
the energy of the model (13. ip . If the coupling constant a is restricted to the set 

-^<ct<0, (3.15) 

then the inequalities (I3.14p yields V 2 iV > and R > 0. The former implies that the flux 
of N at infinity is nonnegative. The positive energy theorem of GR j2H 122] establishes 
that an everywhere positive R gives a positive ADM energy. Therefore, the energy of the 
model, given in (I2.23p . 

E = E ADM - 2a<$> N . (3.16) 
is nonnegative when a is restricted to the range f)3.15p . 



3.2 The IR effective action 

Now we move to the w — > limit to extract the properties of the IR effective action of the 
complete Hofava theory at the A = 1/3 value, which has the potential V^ 2 -* = —R — aaia 1 . 
This effective theory has a consistent structure of constraints by itself, which can be 
seen by taking the w — > limit on the constraints/equations for multipliers of the above 
model. We obtain the ~H and C constraints 

U = —i? j inj ~ VgR + a^g (2V i a i + a^) = , (3.17) 
V 9 

C = ^-p^TTij + IVdNR - (2V iC L l + (2 - a/2)a ia % ) = , (3.18) 

and we notice that the Eqs. (13.7113. 8 p for a and fi maintain their structure of second order, 
linear, elliptic PDEs. By evaluating the Eqs. (12. 26ft . we obtain again the Eqs. (13. 14ft in 
the same form, hence the result about the positiveness of the energy of the IR effective 
action holds. 

Since the number of propagating degrees of freedom is the same of GR, it is interesting 
to elucidate whether GR can be recovered at the lowest energies from this model of 
Hofava theory. The best way to analyze this point is the perturbative analysis. We 
perform perturbations of the Minkowski background in the way 

9ij = f>ij + hj , n lJ =pij , N = 1 + n . (3.19) 

We decompose the perturbative metric into transverse and longitudinal parts, 

K = h 7 + hfoj ~ dijidkk)' 1 ]^ + dihf + fyhf + dijidkk)' 1 ^ , (3.20) 
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where the boundary conditions are used in the definition of (dkk) -1 - In the above we use 
the shorthand dij = didj and so on. We also decompose pij in the same way of (I3.20p . 
In addition, we impose the gauge in which all the longitudinal sector of the metric is 
eliminated, 

dihy^O. (3.21) 

We start the perturbations by analyzing the constraints of the theory at linear order. 
The linear-order momentum constraint W = (12.131) yields 

d iPij = . (3.22) 

Hence the longitudinal sector of p^ is also eliminated. Constraint (12.91) yields p T = 0. This 
leaves us with only the sector pj^ activated. The perturbative Hamiltonian constraint 
"H = (I3.17P and the constraint C = (13.181) at linear order yield respectively 

d ii h T + 2ad ii n = 0, (3.23) 
du^ + Adun = 0. (3.24) 

With the prescribed boundary conditions, and assuming a ^ 2, there are not solutions 
to these equations other than 

h T = n = 0. (3.25) 

Thus, we see that at linear order the transverse scalar mode of the metric and the lapse 
function are switched off completely. Similarly, the linearized equations (I2.27H2.281) yield 
a = fi = 0. 

After all the linear-order constraints are imposed, the unconstrained second-order 
Hamiltonian ( 12.221) takes the form 

H = J d 3 x (pgpg + \d % h]l "dih^j , (3.26) 

which is exactly the Hamiltonian of linearized GR. Thus, we see that the dynamics 
of linearized GR is smoothly and exactly recovered from the complete, nonprojectable 
Hofava theory at the value A = 1/3. There is not any discontinuity in the degrees of 
freedom since the nonperturbative theory has also two degrees of freedom. 

4 Conclusions 

We introduce a theory of gravitation based on the principle of having the foliation- 
preserving diffeomorphisms as gauge symmetry. The theory we consider can be obtained 
from the Hofava theory with the terms of Bias, Pujolas and Sibiryakov by setting the 
coupling constant A equal to 1/3. The terms of Bias, Pujolas and Sibiryakov must be 
included to obtain a renormalizable theory since they are compatible with the gauge 
symmetry. Unlike other models of Hofava theory, in our model the value A = 1/3 
is protected under quantum corrections because of the presence of a constraint in the 
theory. This implies that A is not actually a running coupling constant. 
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Our Hamiltonian analysis shows that the theory has a consistent and closed structure 
of constraints. It also provides elliptic equations for the elimination of the Lagrange 
multipliers associated to the primary second-class constraints. As a remarkable feature 
of the theory, we found that it has two second-class constraints that are absent in others 
models based on Hofava theory. One of them is 7r = 0, which is a primary constraint of 
the theory. These constraints get rid of the extra mode that arises in A 7^ 1/3 models of 
Hofava theory (with or without the terms of Bias, Pujolas and Sibiryakov). We would like 
to stress that these constraints are always present for any choice of the potential, since 
they are a consequence of the universal form of the kinetic term of the Hofava theory 
and the value A = 1/3. As a result, we get that the theory has two physical degrees of 
freedom, as general relativity has. 

We consider a concrete potential with a z = 3 term, a square-Cotton term, and all 
the z = 1 terms compatible with the gauge symmetry. The square-Cotton term can be 
obtained from the principle of detailed balance [I]. For this model we additionally show 
that its energy is positive definite. 

Another outstanding feature of the theory is that the linear-order perturbative ver- 
sion, around Minkowski space-time, of the IR effective action is physically equivalent to 
linearized general relativity. We obtain the IR effective action as the z = 1 truncation 
of the concrete model, since it includes all the second-order terms (including the Bias, 
Pujolas and Sibiryakov term). Thus, we have obtained a clear and consistent way to 
obtain general relativity in the low energy limit of Hofava theory. 

We think that our model is a good candidate for a renormalizable theory of quantum 
gravity. Since the quantization must be performed on the constrained submanifold, the 
second-class constraint it = must be preserved and this avoids A to move from the 
value A = 1/3. However, one should be careful with the treatment of the second-class 
constraints when performing the quantization of the theory. 
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